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Since its discovery, the superconducting phase in water-intercalated sodium cobaltates 
Naa;Co02 • yH2 (x ^ 0.3, y ^ 1.3) has posed fundamental challenges in terms of experimental 
investigation and theoretical understanding. By a combined dynamical mean-field and renormaliza- 
tion group approach, we find an anisotropic chiral d + id wave state as a consequence of multi-orbital 
effects, Fermi surface topology, and magnetic fluctuations. It naturally explains the singlet property 
and close-to-nodal gap features of the superconducting phase as indicated by experiments. 



PACS numbers: 74.20Rp, 74.25.Dw, 74.70.-b 

Introduction. Initiated by the discovery of the 
cuprates, the search for new materials exhibiting uncon- 
ventional superconductivity has become one of the major 
branches in condensed matter physics [1 . A particularly 
exciting idea is the concept of chiral superconductivity, 
where the Cooper pair condensate breaks parity and time 
reversal symmetry and gives rise to interesting edge mode 
phenomena of the bulk-gapped superconductor. As one 
way to accomplish such a scenario, the lattice can act as 
a custodial symmetry to ensure the exact, or, in the pres- 
ence of disorder and nematic fluctuations, approximate 
degeneracy of different superconducting instabilities. In 
such a case, the degeneracy is linked to higher dimen- 
sional irreducible representations of the lattice symme- 
try group, and a chiral superposition of superconducting 
states can be energetically favorable below Tc. Unfor- 
tunately, for the square lattice and its C/^y group, there 
is no such representation for singlet Cooper pairs, which 
in the majority of materials is found to be the generic 
sector for superconductivity. This, however, changes for 
hexagonal systems, where the E2 representation of the 
Cqv lattice symmetry group implies the degeneracy of 
the dx2_y2 and dxy wave state at the instability level [2], 
which can yield a chiral d-\-id singlet superconductor. 

In many respects, the water-intercalated sodium cobal- 
tates NaajCo02 • yH20, with a superconducting dome for 
X ~ 0.3, y ^ 1.3 dit Tc = 4.5K [3 , have been inter- 
preted as the natural generalization of a square lattice 
of copper oxide in the high-Tc cuprates to a triangular 
lattice of cobalt oxide: The electronic structure can be 
assumed effectively two-dimensional due to the intercala- 
tion, and superconductivity emerges as a function of suf- 
ficient Na doping in proximity to magnetic phases. The 
experimental evidence, however, remained inconclusive 
for a significant amount of time, which did not allow to 
draw substantiated conclusions on the nature of the order 
parameter [4]. For example, previous ambiguous indica- 
tions from Knight-shift measurements for polycrystalline 
samples have only later been clarified by single crystal 
measurements [5 , which showed the singlet property of 
the superconducting phase. Similarly, the nodal charac- 



ter of the order parameter has remained a contentious is- 
sue. Early /iSR [6 as well as magnetic penetration depth 
measurements [7 , have shown evidence against a homo- 
geneous gap and have been interpreted in favor of line 
nodes. By contrast, the latest specific heat studies 
advocate a two gap scenario for the cobaltates with one 
comparably small and another slightly larger gap to fit 
the data. This reminds us of similar discussions for the 
iron pnictides, where it is likewise complicated to distin- 
guish a nodal from a strongly anisotropic gap [9 . 

In this Letter, we develop a microscopic theory for 
the nature of superconductivity in the sodium cobaltates 
which is consistent with the experimental findings. In- 
spired by the resemblance to the cuprates, the earliest 
theoretical proposals employed a phenomenological RVB 
theory for the cobaltates [10] supplemented by slave bo- 
son mean- field calculations |TT]. A major challenge from 
the beginning has been the choice of an adequate low 
energy kinetic theory for the problem: While ARPES 
measurements only observe one Fermi pocket centered 
around F in the hexagonal Brillouin zone [12-14 , band 
structure calculations indicate the presence of additional 
Cg pockets [15 . The absence of the pockets in the ex- 
perimental cobaltate scenario has been assigned to sur- 
face effects [16], disorder [17 , and electronic correla- 
tions |T8ti2Q^. It suggests that whatever the microscopic 
theory for superconductivity in the cobaltates may be, it 
should involve a low-energy kinetic theory which explains 
the experimental evidence from a single-pocket scenario. 
Such an effective model has been developed by a com- 
bined dynamical mean-field and cluster approximation 
approach by Bourgeois et al. [211 [22], which is the start- 
ing point of our investigations. Employing multi-orbital 
functional renormalization group (fRG) |9l [23ti26] to ob- 
tain an effective interaction profile for this model, we 
find a rich phase diagram for the sodium cobaltates with 
an anisotropic d + zd-phase in the relevant doping regime. 
The strong anisotropy of the superconducting gap can ex- 
plain the experimental evidence; it follows from the inter- 
play of multi-orbital hybridization, Fermi surface topol- 
ogy, and frustrated magnetic fluctuations in the sodium 
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FIG. 1. (Color online), (a) Effective band structure resulting 
from ^ with t = O.leV, t' = -0.02eV and D = 0.105eV. 
(b) The van Hove singularity is visible in the density of states 
(inset: crystal structure of the C0O2 layers). The Fermi sur- 
face is shown at x = 0.1 (c), x = 0.2 (d), and x = 0.3 (e). The 
Fermi surface colors indicate the dominant orbital weights, 
(d) sketches the division of the Brillouin zone into 96 patches 
used in fRG, (e) depicts the nesting vectors Qn- 
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Cobaltate effective kinetic model Following the work 
by Bourgeois et al. [2TJ|22], all Cos^^ and 02p orbitals are 
taken into account in a finite cluster calculation which is 
then mapped to a three-orbital model and fitted against 
X-ray absorption and ARPES data [22 . The resulting ef- 
fective model is obtained from dynamical mean-field the- 
ory calculations, which adequately take into account the 
self-energy effects at a single-particle level. It exhibits 
strongly hybridized orbitals formed by an effective t2g 
manifold {dxy^dyz^dzx) per site on the triangular Co su- 
per lattice, a finding which is also consistent with ARPES 
polarization measurements [27 . The Hamiltonian reads 
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where fiiaa = Saa^iaa' ^iaa dcuotcs the clcctrou 
creation operator of spin a = t, | in orbital a = 1,2,3 
at site i. t represents the hopping mediated by 02p^, t^ 
corresponds to a direct Co-Co-hopping, D is the crystal- 
field splitting, and /i the chemical potential. We set 
t = O.leV, t' = -0.02eV, and D = O.lOSeV EH |22]. 
The bandwidth of the effective model is ~ 0.6eV. This 
is a factor 3 smaller than LDA calculations predict 
(1.6eV Hg or 2.0eV [28 ). This effective three-band 
model resulting from ([T]) yields one band intersecting 
the Fermi level (Figjl^). A van Hove singularity oc- 
curs at a doping level of x ~ 0.09 (Fig{l]3). The Fermi 
surface contains one hole pocket around F, i.e. the 
center of the Brillouin zone (Fig{l]3). All three hy- 
bridized orbitals contribute to the Fermi surface and 
each has two antipodal dominant regions, indicated by 
red/green/blue dots in Figjlj^-e, corresponding to the 
doping X = {0.1,0.2,0.3}. At x ^ 0.28, the nesting 
of the Fermi surface is optimal, with three nesting wave- 

vectors ^ ^{ ( A O) ^ ^1)' -^)} 

(Fig. [1^). The interaction part of ([T]) features intra- 
orbital Coulomb interaction /7i, inter-orbital Coulomb 
interaction /72, Hund's rule coupling Jh, and pair hop- 



ping Jp. As the bandwidth in this effective model is 
smaller than in bare LDA calculations, the interactions 
strengths from bare ab-initio calculations likewise have 
to be regularized. We set Ui = 0.37eV, U2 = 0.25eV, 
and Jh = Jp = O.OTeV. 

We now proceed by investigating the Fermi surface in- 
stabilities of ([1]) through fRG, where we use the effec- 
tive multi-orbital band structure as the initial starting 
point [9l [23H26] . Through renormalization, we obtain an 
effective low-energy theory of the scattering vertex which 
exhibits superconductivity. The pairing 2-particle vertex 
V^^{k,q) is then decomposed into eigenmode contribu- 
tions which correspond to the different superconducting 
form factors V^^{k,q) = ^icf'^ f{k)f{p). cf^ signals 
the strength of the instability and hence allows to iden- 
tify the superconducting phase adopted by the system. 
We employ multi-orbital temperature-flow fRG [25 to 
take into account the interplay between ferromagnetic 
and antiferromagnetic fluctuations as well as the multi- 
orbital character of the sodium cobaltates. Note, that we 
must avoid the double counting of self-energy effects that 
have already been included to obtain this band structure. 
As such, we intentionally do not take into account self- 
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FIG. 2. (Color online). Phase diagram of the model ([T]) 
as function of doping x and Ui/U2- There are four phases: 
d + id- wave superconductivity {d + id SC, blue), weak ferro- 
magnetism (weak FM, green), /-wave superconductivity (/ 
SC, yellow), and a phase with competitive spin-density wave 
and (i + id- wave superconductivity (SDW / d-\-id SC, purple 
and blue shaded). 



energy effects at the single-particle level which emerge 
during the RG flow. In total, our procedure is still an 
approximation, in the sense that we first renormalize the 
single-particle level via a DMFT approach and then sep- 
arately investigate the renormalization of the scattering 
vertex through fRG which gives rise to Fermi surface in- 
stabilities. This, however, is justified because the scat- 
tering vertex evolution under RG is only significant in 
the immediate vicinity of the Fermi level. 

Phase diagram. From an itinerant viewpoint of Fermi 
surface instabilities, an important feature of ([T]) is the dif- 
ferent doping location of the van Hove singularity {x = 
0.09) and the optimally nested Fermi surface {x = 0.28). 
Note that both locations are coincident and much less re- 
vealing for a triangular lattice tight-binding model with 
only nearest neighbor hopping [29 . Accordingly, for 
small doping, the phases are determined by the large 
density of states at the Fermi level combined with rather 
weak nesting corresponding to finite momentum transfer, 
which in total promotes dominant zero momentum parti- 
cle hole scattering (labeled weak ferromagnetism (FM) in 
Fig.[2|. With increased doping, it depends on the ratio of 
U1/U2 whether the system favors triplet /-wave or singlet 
d + id-waye superconductivity. For the former, the sys- 
tem exhibits nodes along the Fermi surface and follows 
the gap function f {k) = sm{ky) - 2 cos(^^) sin(^) 
(Fig. |3^,b). /-wave is preferred for enhanced Ui/U2^ 
as it is seeded by spin alignment stemming from ferro- 
magnetic fluctuations, which are reduced by • In the 
case of preferred d-wave superconductivity, we find two 
degenerate instabilities associated with the form factors 
depicted in Fig. [Sfc, which relate to the leading harmon- 



ics dx2_y2 (k) = 2cos(kx) — cos(- 



) — cos( 



and d^y ik) = cos(^^±^) - cos(^^^^). (Note that 
throughout parameter space, we always find higher har- 
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FIG. 3. (Color online). Superconducting form factors from 
fRG mean-field decoupling of pairing channels (solid lines) 
compared to analytic leading harmonic solutions (dashed 
lines): (a,b) /-wave, (c,d) d + id- wave. The representa- 
tive parameters are (a,b) U\/U2 = 1.4, x = 0.18 and (c,d) 
U1/U2 = 1.0, X = 0.14. (b) and (d) show a color plot of the 
gap size Ao along the Fermi surface (Eq. |2|. 



monic contributions in the (i-wave sector [30^ to be irrele- 
vant.) The system could generically form any linear com- 
bination di + 6*^(^2 of both d-wave solutions which must 
be degenerate at the instability level as protected by lat- 
tice symmetry. A mean-field decoupling in the SC pairing 
channel and minimization of the free energy as a func- 
tion of the superposition parameter, can be rephrased by 
satisfying the self-consistent gap equation [2 



2E{k) 
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We always find d -\- id to be the energetically preferred 
combination. This is rather generic in a situation of de- 
generate nodal SC order parameters, since such a combi- 
nation allows the system to avoid nodes in the gap func- 
tion (Fig.jsji) and maximizes condensation energy. Note, 
however, that the relative energy gain between di d -\- id 
state and a different possible solution such as nodal single 
dx2_y2 varies significantly depending on the microscopic 
setup. For example, the condensation energy gain from 
d-\-id for lower doping in Fig. ^ as compared to a single 
d wave solution will be higher than for larger doping in 
Fig. [if. 

Anisotropic regime at x ^ 0.3. The effective model 
in ^ is quantitatively most accurate in the doping 
regime of the superconducting dome of the cobaltates. 
There, we find d -\- id superconductivity with strong 
SDW fluctuation background (Fig. [2|. The enhance- 




FIG. 4. (Color online) Change of gap anisotropy in the d-\-id 
superconducting phase as a function of doping, (a-c) x = 0.14. 
The dominant density of states is strongly peaked at the 
warped edges of the Fermi surface and singles out specific 
particle hole scattering channels which result in a homoge- 
neous gap Aq. (d-f) X — 0.28. Various scattering channels are 
enhanced due to reduced Fermi surface warping and more ho- 
mogenous distribution of density of states, yielding a strongly 
anisotropic gap. 



merit of magnetic fluctuations is also observed in experi- 
ment, anticipating the metal insulator transition regime 
at X ~ 0.5. From the viewpoint of Fermi surface topology, 
this is due to improved nesting conditions for a larger part 
of Fermi level density of states (Fig. [l^). As the density 
is also more homogeneously distributed along the par- 
allel sides of the hexagonal Fermi surface, however, this 
leads to an increased bandwidth of enhanced particle- 
hole channels in the RG flow, which eventually yields an 
enhanced anisotropy in the seeded pairing channel which 
gives rise to d -\- id superconductivity. This trend for 
increasing doping x is illustrated in Fig. [4] At lower dop- 
ing, the more warped Fermi surface clearly singles out 
the particle hole channel according to scattering between 
the Fermi surface edges, which also possess the dominant 
fraction of density of states at the Fermi level (Fig.[4^-c). 
As doping is increased, the reduced warping allows the 
effectively one-dimensional parallel sections of the Fermi 
surface to drive more particle-hole channels. The result is 
an evolution of the seeded superconducting phase from 
a homogeneous to an extremely anisotropic d -\- id gap 
(Fig#-f). 

In Fig. [Sj we have plotted the degree of gap anisotropy 
in the d-\-id phase for the same range of parameters as for 
the phase diagram in Fig. [2] We use the variance of the 
gap function, divided by the mean: ri = Notably, 
the physically relevant regime for the cobaltates coincides 
with the strongest gap anisotropy, we observe. As seen 
in the inset of Fig. [5) this yields a gap structure along a 
single Fermi surface which can be well characterized by 
a very small gap along the sides and a large gap at the 
edges of the Fermi surface. It can explain the two gap 
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FIG. 5. (Color online). Gap anisotropy r] = ^^e 
d + id- wave phase (axis annotations as in Fig. [2]). Blue (red) 
regions indicate a rather homogeneous (anisotropic) d-waxe 
gap. The inset shows the plot of the gap function on the 
Fermi surface at physically sensible doping: Ao exhibits close- 
to-nodal (blue) and larger (red) gap regions along the Fermi 
surface. 



scenario from specific heat measurements [8^ , where both 
effective gap scales originate from a single pocket. 

Conclusion and perspectives. We have shown that the 
experimental evidence of a close-to-nodal singlet super- 
conducting state in the sodium cobaltates can be de- 
veloped from a microscopic model taking into account 
the multi-orbital nature of the electronic scenario. The 
anisotropic d -\- id superconductor we find is a combined 
effect of magnetic fluctuations, specific Fermi surface 
topology at the corresponding Na doping, and multi- 
orbital effects. To our conviction, this constitutes the 
sodium cobaltates to be one of the most promising can- 
didates for a chiral singlet superconductor, to be fur- 
ther studied experimentally and theoretically. The Tc 
might allow for laser-ARPES studies in the supercon- 
ducting phase, along with more careful investigations of 
time-reversal symmetry breaking than it has been pur- 
sued by now. Likewise, the role of lattice distortions and 
disorder can be interesting to consider, as there might 
be a transition from d -\- id to d when the custodial Cq^ 
symmetry is sufficiently broken. Finally, the desire for 
an ideal chiral singlet superconductor might also warrant 
further work to optimise its two-dimensional character, 
as has recently been reported for intercalated iron-based 
superconductors [31] 
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